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1. Introduction 

In the chiral limit of QCD, obtained by setting the light quark masses to zero, the global sym¬ 
metry of QCD is chiral SU (2) x SU (2). This is homomorphic to a 0(4) symmetry. The universality 
conjecture then leads us expect that the critical indices and critical amplitude ratios in chiral QCD 
should be the same as for 0(4) magnetic systems [|^. 

In recent times, the chiral condensate of QCD (which is equivalent to the spontaneous magne¬ 
tization of the magnet), its derivatives with respect to the quark mass (equivalent to the magnetic 
susceptibility), and their scaling towards the chiral limit, have been studied with somewhat am¬ 
biguous results 0. 

The scaling of the internal energy and the specific heat of 0(N) symmetric systems is intricate, 
as is known from the phenomenology of liquid Helium. However, with the extensive lattice QCD 
computations now available on quark number susceptibilities in QCD at zero baryon density |^], 
their important role in heavy-ion collisions 0], and their connection with temperature derivatives 
of the free energy, it is important to initiate the scaling analysis of these quantities. That is the 
purpose of this talk. 


2. Scaling and the limits of universality 


The thermodynamics of QCD is characterized by a free energy, which is a function of some 
number of intensive control parameters. These could include the temperature, T, and the baryon 
chemical potential, jj.. If the pion mass were exactly vanishing, then QCD would have a 0(4) global 
chiral symmetry. Since we are interested in real QCD, where the pion is not massless, an explicit 
chiral symmetry breaking parameter is needed. This is the quark mass, m, which plays the same 
role in QCD as a magnetic field does for the 0(4) magnet. 

Near the critical point one can decompose the free energy, F{T,m), into the sum of two terms. 
One of these is a regular part, Fr{T,m), and the other is a singular part, Fs{T,m). Fr{T,m) is Taylor 
expandable around the critical point, T = Tc and m = 0 with some large radius of convergence. 
The modern theory of critical phenomena starts from the observation that the most singular part is 
a scaling function 


Fs{T,m)=t^ “<I>(t), where 



t 


(2.1) 


where we have chosen the scaling variables t and T to be dimensionless, Tc is the critical temper¬ 
ature, Mo is any mass scale which remains finite in the chiral limit, and a, A = j35 are critical 
exponents. The function defined so, is universal, in the sense that whether we examine 

an 0(4) Heisenberg magnet, QCD, a non-linear sigma model of pions, or the Nambu-Jona-Lasinio 
(NIL) model, the Fs{t, z) we obtain from all of them are the same. As a result, the universal prop¬ 
erties of thermal QCD know nothing about QCD, aside from its chiral symmetry. However, the 
various models differ in Fr{T,m), so this is the piece which gives information about the actual 
degrees of freedom involved in the QCD phase transition. 

Even if the magnitude of Fs is comparable to iy, since it is singular, its effect may be enhanced 
by taking sufficient number of derivatives. For example, the specific heaf, Cy and hence 
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Model 

Example 

i8 

5 

a 

Ref 

O(oo) 


1/2 

5 

-1 

§] 

0(4) 

chiral QCD 

0.380 

4.86 

-0.2268 


0(3) 

? 

0.365 

4.79 

-0.115 

0 

0(2) 

liquid He 

0.349 

4.78 

-0.0172 

0 

0(1) 

liquid-gas 

0.325 

4.8 

0.11 

0 

MFT 

1/2 

3 

0 



Table 1: Critical exponents of 0(N) models in three spatial dimensions. The exponent a is obtained from 
the other reported exponents using scaling identities. 0(1) should be taken to mean the Ising model. MFT 
stands for mean field theory. There are no known examples of 0(3) models, since real ferromagnets have 
relevant terms which break this symmetry. 


diverges at T = Tc, provided that a > 0. As one can see from Table H, this is true of the Ising 
model. However, for all other 0{N) models a < 0, and, as a result, the singular contribution to the 
specific heat exactly vanishes for T = Tc. 

This seems to contradict our knowledge of the specific heat of liquid He, which is in the 0(2) 
universality class and has a cusp at the critical point. The resolution of this puzzle comes from 
noticing that the peak of Cy is only finite, and hence is regular. It is the shape which is singular. 
So the specific heat has to arise through a playoff between the singular and regular parts. In fact, 
a very precise microgravity experiment has been done over the range — Tc\<2 nK |§], and the 
results fitted to the formula 

Cy=Ar + r^{B + Ct-^), ( 2 . 2 ) 

where comes only from /y. A' is a possible correction-to-scaling exponent. Ar is positive and 
B is constrained to be negative. It can be shown that B can be negative without violating the 
thermodynamic consistency criterion that Cy > 0. A result of the microgravity experiment is that 
a = —0.01285(38). Interestingly this is in disagreement with the careful work of 0. 

This mechanism also works for QCD and other 0(N) symmetric models. For these we may 
write 

Cy{T,m)=Ar + ^—^{t,T). (2.3) 

Tc and Ar are non-universal, and change from QCD to various effective theories for it, but the 
exponent a, and the regular function 'T (which may be written in terms of the scaling function <I>, 
and its derivatives, <!>' and <I>", if desired) are universal. So the shape of the specific heat cusp is 
universal but its height and width must be determined in QCD. Furthermore, these two parameters 
are good tests of possible effective models, since a bad model of QCD will not reproduce its non- 
universal properties. 

An interesting statement about the scaling of Cy with mass arises from this. Suppose we 
succeeded in measuring (on the lattice) Cy for QCD with various different light quark masses. By 
plotting the data as a function of t and scaling Cy appropriately, can we observe scaling in the form 
of data collapsing on to an universal curve? Clearly, there are no singularities of the free energy if 
T is varied around Tc at fixed non-vanishing m. As a result, taking t —> 0 and T ^ 0 simultaneously 
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Figure 1: Data collapse obtained in an MFT treatment of the NJL model when keeping only the data for 
T > 50. Tiny violations of scaling are visible; these can be controlled by increasing the cut on T. The value 
of this cut depends on the choice of Mq (here it was taken to be Tc), and the renormalization scheme, when 
going beyond MFT. Also, since this value is not universal, it could be different in QCD. Only data for T <Tc 
is used in this plot to avoid having to subtract a large regular part, as discussed in the text. 


will not reveal scaling. Instead, one must take the limit m ^0 first and T ^ Tc next, which means 
that one must take T —)• oo first and then t —)■ 0 in order to see data collapse. 


One can test this in the NJL model even at tree level, i.e., in the MFT approximation. The high 
temperature limit of this model contains weakly interacting quarks, so the regular contribution to 
Cy actually increases fairly rapidly with temperature. As a result, one may miss the pseudo-critical 
behaviour in Cy unless the temperature range is scanned finely to discover a peak sitting over a 
rising background, or the free quark contribution is subtracted to make the peak stand out over a 
falling background. This difficulty would also occur in QCD [IF], but not in the 0(4) Heisenberg 
magnet. 

With this MFT one sees that data collapse is possible when one plots Cy against t provided that 
one selects only x > 50. This is sufficient to ensure that for any finite m one does not approach 
t = 0 too closely. Figure || shows that one may relax the condition T —> oo provided one is willing to 
tolerate small enough violation of scaling. Since experimental data or Monte Carlo computations 
come with errors, it should be possible to tune the cutoff on T in order to find the scaling curve 
within the errors. 


We end this section with a remark about the scaling fields. In making use of effective models 
to study universal properties of QCD, most works make the assumption that the scaling fields of 
the effective theory {T and m) are identical to those of QCD. Whether or not this assumption is 
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Figure 2: Data collapse obtained for Xb using the results of 0^ when keeping only the data for T > 0.04 
with 0(4) exponents. The data is plotted as a function of \ — T/Tc instead of t because the regular parts on 
the two sides of Tc are have not been removed. The colour coding corresponds to the values of irl^lm^p given 
in the figure legend. 


correct can be tested, but, to the best of our knowledge, such tests have not been performed. 


3. Relevance to the phase diagram of chiral QCD 


The phase diagram of chiral QCD can be extended to finite baryon chemical potential, jx. 
Since this scaling field preserves fhe 0(4) symmefry, fhe critical poinf of chiral QCD af /r = 0 gels 
slrelched info a line. The global CP symmefry of QCD implies fhal F{T,m = 0,IJ.) = F{T,m = 
0, —ft), so 

r,(ft) = r, + itcft2 + .... (3.1) 

The curvalure has been sludied in laffice QCD for abouf a decade. Differenl deferminafions agree 
roughly on ifs value iflTl]. 

If one assumes fhal ft enters the scaling function of eq. ( |2.1[ ) only through the dependence of 
Tc on ft as given in eq. (0), then, as in [|I^, we can write 



fi=0 


Tc 

TKdjX^ 




p=0 


(3.2) 


Qualitative evidence for this relation between derivatives was obtained very long back [[T^. The 
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relation above implies a connection between Cy and the fourth order baryon number susceptibility- 


xtiTdn) = 


d‘^P{T,m,n) 




11=0 




Cy 


(3.3) 


The last expression for comes from retaining only the most singular contribution. This suggests 
a scaling test of Xb similar to that for Cy. 


This test can be performed with the results of 114]. For the treatment of lattice QCD computa¬ 


tions we may replace ui/Mq in the definition of T by so that t, T and Xb renormalized 


t^lrn^p, so that t, T and Xb '• 
quantities. In Figure § we plot the full measured TcXb/T■ Data collapse should be expected in the 
region where the contribution from the singular part dominates. 

The regular parts at temperatures well above and well below Tc are expected to be different, 
since effective theories in these two regions are the hadron gas model {T <C Tc) and the weak 
coupling expansion of QCD {T S> Tc). We have plotted Figure § to show these two branches 
separately. It seems that in the region t <0.1, the differences in the regular parts may be neglected 
within the precision of the data. One particular implication is that the gas model should not work 
for Xb within 10% of Tc. 

While this gives us a first test (in this sector) of scaling at T = 0 at surprisingly large quark 
masses, the current errors are large. Improvement in errors would allow us to test scaling better, and 
also to test the importance of the variation of the regular part of these quantities with approximately 
10% change in t. These requirements set benchmarks for future measurements of Xb- 

It is clear from eq. ( ^T| ) and the values of a in Table |I] that derivatives of Cy with respect to T 
would diverge in the vicinity of the critical point. As a result, one should be able to observe scaling 
of the higher order baryon number susceptibilities. For example, the sixth order quantity, Xb would 
be universal, and is likely to have a shape similar to that shown in Figure |3[ It would be interesting 
to test this in future, when improvement in statistics makes these tests significant. 


4. Some remarks 

It is interesting to recall that before the modem theory of critical phenomena was developed, 
the Ehrenfest classification of phase transitions was in vogue. This attempted to define orders of 
phase transitions according to which derivative of the free energy diverged. In the case of 0(N) 
models one sees very clearly that such a classification runs into trouble. On examining the chiral 
(magnetic) susceptibilities, one comes to the conclusion that the QCD transition is of second order. 
However, on examining Cy one comes to the conclusion that the same transition is of third order, 
since Cy does not diverge, but its derivative with respect to T does. We realize today that the 
differences are due only to the value of associated critical index. 

We conclude by reiterating the importance of scaling tests such as that suggested here. They 
constitute a new domain of tests of the universality hypothesis in the context of QCD. Not only 
is this important in its own right, but also serves to put bounds on the region of applicability of 
models such as the hadron gas model. This model is a mixture of ideal gases and hence contains 
no singular part, whereas 0(4) universality is based entirely on the singularity due to pions in the 
chiral limit. Since these are mutually exclusive descriptions of the free energy, the success of one 
rules out the other. 
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Figure 3: In the chiral limit Xb diverges with critical exponent 1 + a and is proportional to the temperature 
derivative of Cy At finite rn^jm^ these divergences would be rounded off as shown. 


We thank Deepak Dhar for discussions. 
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